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We investigate the collective excitations of a low temperature dilute gas mixture that consists 
. of a Bose-Einstein condensate and a Fermi-gas that is a normal (i.e. non-superfluid) Fermi-liquid. 

, We find that the BEC-mediated fermion-fermion interactions, as a consequence of retardation, can 

become repulsive and support a zero-sound mode that is essentially undamped. In addition, we find a 
damped zero-sound mode that can be described as a BEC-sound mode modified by fermion mediated 
boson-boson interactions, and we derive its decay-rate caused by Landau damping. We study the 
. mode structure of these excitations and find avoided crossing behavior as well as a termination point. 

O ' The collective mode dynamics also reveals that phase separation sets in when the fermion-mediated 

a CD ■ 
boson-boson interaction destroys the stability of the homogeneous BEC. We estimate the time and 
length scales of the onset of the phase separation, and we discuss the feasibility of experimentally 
probing these consequences of mediated interactions. 

C/5 



O 

o 



> 
o 



I 



PACS numbers: 05.30.Jp, 03.75.Kk, 32.80.Pj, 67.90.-l-z 



I. INTRODUCTION 



. Current atom trap experiments create quantum degenerate Fermi-gas systems 0, and Fermi-Bose gas mixtures 
' in sympathetic cooling experiments 0. The sympathetic cooUng of fermion atoms utihzes bosonic atoms at lower 
temperature, often condensed into a Bose-Einstein condensate (BEC). At the time of writing, the fermions do not 
quite enter the Fermi-liquid regime before the BEC is separated out and the fermions arc further cooled evaporatively 
[3|; however, unambiguous Fermi-hquid/BEC mixtures are attainable Fermi- liquid behavior sets in when the 
energy can be expanded up to second order in the (quasi-) particle occupation numbers, the heat capacity varies 
hnearly with temperature, and the Pauli exclusion principle sufSciently inhibits incoherent fermion scattering to 
validate the assumption of coUisionless dynamics. A dilute two-component Fermi-gas enters this regime when the 
^vq ' temperature drops below 10% of the Fermi temperature |^. Strong interaction effects further lower the Fermi- liquid 
QQ range: experiments reported an upper bound of 2% of the Fermi-temperature in the condensed '^He liquid. 

With the creation of cold atom Fermi-Bose quantum liquid mixtures, atom traps can access novel regimes of low 
•/^ , temperature physics. At present, table-top explorations of Fermi-Bose quantum-liquid mixtures are confined to the 
■ condensed '^He-^He mixtures. These experiments revealed an intricate phase diagram, even if the ^He component did 
not reach superfluidity. The critical temperature for ^He superfluidity drops in the '^He 'rich' regime, and depending 
on the temperature and fraction of the '^He isotope, the mixture spontaneously breaks translational symmetry by 
forming ''^He only' regions. The microscopic understanding of this phase separation is, however, complicated due to 
^ the strong interaction which masks or, at least, competes with the quantum statistical effects @. 
r"| ■ General arguments based on Fermi- liquid theory indicated the importance of mediated interactions they are 
Q , the modification of the inter-particle interactions in one liquid due to the presence of the second liquid and directly 
O ■ causes the phase separation by destroying the stability of the homogeneous state. In the static limit, applicable to 
the description of the helium mixtures in which the sound velocity c of the boson superfluid greatly exceeds the 
, , • Fermi velocity vf, the boson- mediated fermion-fermion interactions are attractive and can Cooper-pair the fermions 
rS ' [sl 1^ ITo|. From the experimental measurements in helium mixtures, it appears that the strong interactions greatly 
weaken the mediated interactions f^. 

In contrast, the phase separation of cold atom Fermi- mixtures [TTI | would be amenable to a transparent first- 
principles description [T^ and the boson-mediated fermion-fermion interactions can significantly increase the critical 
temperature for Cooper-pairing |l3j| . In addition, the relative densities of the cold atom mixtures can reverse the 
Fermi- velocity vf to Bose sound velocity c ratio to yield c/fp < 1, which is not attainable in the helium mixtures. 
Under those conditions, satisfied in most sympathetic cooling experiments, the static description breaks down and we 
show that retardation radically alters the nature of the boson-mediated fermion-fermion interactions. 

There are other aspects that motivate cold atom Fermi-liquid BEC studies. For example, the range of the BEC- 
mediated fermion-fermion interactions, which is the BEC-healing length, can exceed the inter-particle distances. Cold 
atom experiments with fermion-BEC mixtures can then go beyond the usual cold atom paradigm of short-range. 
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contact-like inter-particle potentials while examining similarities with other Fermi-liquids. In particular, the form of 
the static BEC-mediated interaction suggests an intriguing connection between the anticipated cold atom Cooper- 
paired Fermi-Bose mixtures and the BEC-exciton superfluid that has been pursued in condensed matter for more 
than two decades The BEC-mediated interaction potential is of the Yukawa form which resembles the 

electron-hole Coulomb attraction of an exciton at distances smaller than the Yukawa range. 

The realization of s-wave Cooper-pairing in fermion-BEC mixtures, will, however, require much experimental skill. 
This pairing scheme requires the simultaneous trapping of fermion atoms in two distinct hyperfine states at nearly 
equal densities, a considerable challenge in the presence of a BEC: A dependence of the boson-fermion interactions on 
the fermion spin introduces a difference in the local spin densities even if the total number of fermion atoms is evenly 
divided over the two spin states and both spins experience identical trapping potentials. 

In this paper, we focus on a dilute mixture of BEC and one-component Fermi gas that remains normal. We 
present a theoretical study of the small amplitude oscillations around a spatially uniform ground state. First, we 
show that the resonant dynamics of the mediated interactions introduces novel features. If c < vp, the mediated 
interaction supports a long-lived essentially undamped zero sound mode as in a two-component Fermi-liquid case, 
which is incompatible with the effective attraction predicted by the static treatment. On the other hand, if vp < c, 
the mixture supports an additional mode that is the remnant of the BEC-sound mode modified by a cloud of fermion 
excitations that introduce fermion mediated boson-boson interactions. This second mode is damped. Both oscillations 
have been discussed by Yip '16] from a calculation of man y-b ody response functions. Several calculations of zero- 
sound modes in trapped fermions have been also reported in \IT\ . We analyze the homogeneous BEC-fermion mixture 
using the semi-classical description, and obtain a single self-consistent Fermi-liquid-like collective mode equation. The 
semiclassical description in terms of the spatially oscillating deformation of the local Fermi-surface elucidates the 
zero-sound nature of the fermion oscillations in both modes, revealing that the dominant decay mechanism of the 
BEC-like sound mode at long wavelengths is Landau damping. 

The dispersion relations of the two collective modes, the undamped and damped modes, avoid crossing each other 
in the wavenumber interval where the undamped mode gradually changes its character from a pure Fermi-liquid zero 
sound mode to a Bogoliubov excitation of the BEC. The damped mode terminates at a certain wavelength. The 
damping rate at the termination point remains significantly smaller than the excitation energy so that the mode 
represents a good quasi-particlc for all wavenumbers and the termination point should be observable. 

The phase separation is manifested by the collective mode dynamics. The instability corresponds to collective mode 
eigenvalues of purely imaginary value. Our analysis reveals the mechanism of the instability, and gives the length and 
time scales on which the phase separation sets in. 

First, we discuss the mediated interaction of the fermion-BEC in a static limit in Sec. II. Then, the dynamical 
effects are obtained in Sec. Ill, followed by the mode analysis in Sec. IV. Finally, experimental aspects are examined 
in Sec. V. This is an extension of the work previously presented in Ref. 18]. 

II. MEDIATED INTERACTIONS IN THE STATIC LIMIT {lo 0) 

In most mixtures, the mediated interaction shifts the frequency of the existing mode. In the single-component 
fermion BEC mixture, the BEC-mediated interactions actually support collective Fermi-liquid like oscillations. Inter- 
actions amongst identical Fermions are prevented by Pauli-exclusion, so that any fermion-fermion interactions in the 
single component fermion BEC mixture are necessarily BEC mediated. 

In this section we illustrate BEC-mediated interactions by considering the simplest system to exhibit this phe- 
nomenon: two indistinguishable, stationary fermions of mass mp are embedded in a dilute BEC of the particle 
density and particle mass uib in equilibrium located at position xi, and X2 . The mediated interaction obtained 
in a static system highlights the dynamical effects that we discuss in Sec. IIIII Without the impurity atoms, the BEC 
would be homogeneous. The bosons interact by short range interactions, which we describe by a pseudo-potential 
XbbS{t — r') with interaction strength. Abb = 47r/i^aBB/'7iB, that is proportional to the boson-boson scattering length 
Abb- Similarly, we assume that the interaction between the fermion and BEC-boson atoms is well described by 
X]j -^Bpf^lxj — r), with an interaction strength proportional to the impurity-boson scattering length obf and inversely 

proportional to the effective mass of the fermion-boson system ttibf, where m^p = rrip^ + , Abf = 27rfi^aBF / "t-bf- 
The pseudo-potential description presupposes that the scattering lengths are shorter than the average inter-particle 

distances, aBB,aBF ^ {pb) 

This system lends itself to a particularly transparent description of the mediated interaction. The BEC-mediated 
interaction is the modification of the mean-field energy experienced by the fermion at Xi due to the other fermion at 
position X2. The energy shift is caused by the change in BEC-density /5b(xi) due to the interaction of the other atom 
with the BEC. Within the mean- field description, the equilibrium BEC field in the presence of the X2 atom follows 
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from the time-independent Gross-Pitaevskii equation 



MB'?^(r) 



2m} 



'- + \bbW)\ 



(r) + \bf5{v - X2)0(r) , (1) 



where 0(r) represents the condensate field (j) — ("0) with ps(r) = I^Cr)]^ and /ie — ^bbPb is the chemical potential of 
the BEG ensuring that the density tends to at large distances from X2. Since the BEG is in equilibrium, the field 
(j) can be chosen to be real-valued. 

Linearizing the condensate wavefunction by setting w VPb + ^4>^ Eq.(^ becomes 

[v2-r'lj0(r)--S<5(r-X2), (2) 



where ^ is the healing length of the BEG (also referred to as the coherence length) , ^ = 1/ IQtt p'^a-g.B and ^ is the 
strength of the source term, S = — 47r(l -f mB/mp)a-pB\/l^- The solution to Eq. ^ is proportional to the Green 
function of the modified Helmholtz equation, 

1 + -e « . (3) 

ruY J |X2 - r| 

The condition, |(S(/)(r) = C <C 1 suggests that the Yukawa profile is accurate if |xi — X2I > Vg, where Vg = 

aBB(l + "^b/"t^f)/C7 and we have assumed that ^ > Ts- The BEG-mediated fermion-fermion interaction potential in 
the static limit is 



Vr"" (|xi - X2I) = 2Afb Vp^<5</)(xi) 



= _2 f 1 + Afbp^fb^^^^^^^V^. (4) 

V "^F/ IX1-X2I 

The result is an attractive Yukawa potential, familiar from the description of effective interactions mediated by 
massive scalar fields in relativistic field theory For BEG's this result was obtained as well using perturbation 

theory [Hill El. 

From Eq. we find that the average value of the static mediated interaction for a fermion particle in a homogeneous 
BEG-fermion mixture is only a fraction of the mean-field fermion-boson interaction energy experienced by a fermion, 

'^fbPb' ■^ith tl^s fraction being of the order of 2 (1 -I- mB/mp) {p%)^^^ o-fb exp (— However, in the long 
wavelength limit, the Fourier transformed potential, V^°'^(fc), is 

V^^\k)^-^ (5) 

Abb 1 + (fcO 

and comparable to or possibly larger than that of the effective contact potential, Afb, linife-»o V-^'^{k) — — A|b/Abb- 
We note that Eq. © agrees with the long wavelength limit derived in Ref. jlBj using the Landau Fermi-liquid approach 
(which is also valid in describing the strongly interacting systems). 



P,p'->o dps op- 



where Ep is the fermion quasi-particle dispersion, and the derivatives are taken while keeping the occupation numbers 
constant. For the weak-interaction, Ep k, AfbPb +?'^/2™f, and we recover the long-wavelength limit of Eq. Q. 

Similarly, we can determine the fcrmion-mediated interactions in the static limit by calculating the density variation 
caused by a single impurity atom located at x, interacting with the fermions via a short-range interactions. The 
resulting density variation oscillates on the length scale of the average fermion-fermion distance, yielding a density 
profile PF(r) with features known as Friedel-oscillations. At large distances, the density variation (5pF(r) of a three- 
dimensional weakly interacting fermion system oscillates as 

cos(fcF|x-r| + k) 

(5pF r cx 13 7 

|x — r]-^ 

where k-p is the Fermi momentum, ^f = (67r^pp)^/'^, and k is the phase shift. The resulting mediated interaction is 
known as the RKKY-interaction. 
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III. COLLECTIVE MODE DYNAMICS 



While the static limit of the mediated interactions presents a picture of appealing simplicity, most sympathetic 
cooling experiments operate outside the boundary of its validity (which amounts to v-p <C c). In this section, we 
investigate the effects of the dynamics of the BEC-response in collective modes. The BEC density fluctuation propa- 
gates at a finite velocity and the consequent retardation of its response to a fermion density fluctuation can alter the 
physics of the collective oscillations radically as we show below. 

We assume that the boson-boson interactions are repulsive, aes > 0, and the temperature is sufhciently low so that 
almost all boson particles are condensed into the BEC, thus the fermion dynamics can be described as collisionless. 
These assumptions will lead to a description of the dynamics that is effectively temperature independent even though 
it is tacitly understood that the temperature exceeds the critical temperature for fermion pairing. 



A. Fermion dynamics 

The inhibition of incoherent fermion-fermion scattering by Pauli-exclusion principle does not eliminate the effects 
of fermion-fermion interactions. As shown by Landau |20j, collective modes can be driven by the mean- field varia- 
tions that accompany the propagating density and current variations. Zero-sound, which corresponds to a periodic 
deformation of the local Fermi-surface, is the quintessential example of such motion. 

Let us define the single particle fermion distribution function as n(r,p;i), where r is the location and p is the 
momentum. The time evolution of this function is governed by the collisionless transport equation 

9n (r,p,t) dn (r, p, t) , p . ^ ,^ dn (r,p,t) _ 

dt ^^^^ d-r +F(r,P,i)- -0, (8) 

where we have approximated the velocity of the dilute fermion gas as v ~ p/m-p, and F is the gradient of the 

BEC-induced mean-field interaction energy, F = —V [AfbPb (r, t)], and, in Cartesian coordinates, V = . 

Then, we linearize n assuming that n(r, p; t) fluctuates around the equilibrium distribution (p) and the fluctuation, 
which is confined in momentum space to a thin shell in the vicinity of the equilibrium Fermi-surface, propagates in 
real space as a plane wave with wavevector k. We write n (r, p, t) as 

n(r,p,i) =nO(p)+WpPF5(|p|-PF)e'(*'""-"*', (9) 

where pp is the equilibrium Fermi momentum and Up denotes the amplitude of deformation. Similarly, assuming that 
the BEC fluctuation also propagates as a plane wave, /C_B(r, t) = p'^+SpB (r, t) with 6pB (r, t) — Tj^^p'^ exp [i (k • r — ujt)], 
we obtain from Eq. © with Eq.@ that 

k • P k • P r, 

-lujup + i w -I — —AfbPb^Ic, (10) 

In the plane wave propagation of the BEC fluctuation, the complex phase of the BEC-field varies as 56b = 
9k exp [i (k ■ r — tut)] . The assumption of small amplitude oscillations implies that the fluctuation amplitudes 0k 
and 77k ^ 1- Below, we describe their time evolution. 



B. BEC dynamics 



In general, the BEC dynamics is conveniently described in terms of the real valued density and phase fields, pB(r; t) 
and 0B(r, t), the fluctuations of which we introduced above [2]J. These quantities relate to the complex valued 

condensate field as (t>{v,t) — ^■i/'(r, i)^ = ^J'pBiv^e^^'^^'^'*\ The time-dependent Gross-Pitaevskii equation (Eq. 

with /iB — > —ihd/dt) leads to 



i2mB 



= 0. 



With the density and phase fields, Eq. 1)11(1 takes the form of a classical continuity equation, 

dpB{r;t) 



dt 



-V- [pB(r;t)v, 



(11) 



(12) 
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where Vg = {h/m-o) V6'B(r;i) is a superfluid velocity. The real part of the Gross-Pitaevskii equation gives 

96*3 (r,i) -h? —2 / ^-T 

— ii — - = / , V Vpb (r, t) 

dt 2mBv/pB (r,i) 

+ AbbPb (r,t) + Afb/Of (r, i) + imBV^(r;i), (13) 

2 

which is similar to the Landau superfluid velocity equation with the Bohm-potential (oc V ^/pb) added to the usual 
mean-field contributions AbbPb + AfbPf and the superfluid kinetic energy, mBVg/2. 

Linearizing the BEC-fields as in Sec. IIII Al then substituting them into Eqs. (|12|l and H13(l give 

-iw?7k = (14) 
tob 

ihujO-i^ = + AbbPb j ?7k + ^fbSpfm , (15) 

where (5pF,k is the Fourier transform of the fermion density fluctuation that corresponds to Eq.©. The fluctuation 
amplitude of the boson density is then proportional to the fermion density fluctuation and can be written as 

^ ApB 5pvm/Pb 
"^■^ Abb [io/kcY - [1 + [kif] ' ^ ^ 

1 /2 

where c is the sound velocity of the pure BEG, c — (Abb/Ob/'tib) • Note that the sign of the ratio of the boson to 
fermion density fluctuations depends on the frequency w. 



C. Dispersion relation of the collective modes 



Returning to the Fermi- liquid dynamics of Sec. IIII 51 we write Eq. (|10|) as 



Up - 77k cos {ip) 



A 



fbPb 



2eF 



= 0, 



(17) 



where Lp is the angle between the k and p -vectors, and where ep represents the equilibrium Fermi-energy. 

From the definition of n, which implies pp{v,t) — J dPp n(r,p,i), and with Eq. @, we find that the Fermi- 
density fluctuation SpF,k in Eq. (|16|l is proportional to the angular average of the deformation amplitude, 



where (u) = {iny^ J 



Combining Eqs. (|16|l and (|17l) . we obtain 

F 



kvF 



COS {if) 



COS (ip) 



[io/kcY - [1 + {kS,f 



(18) 



(19) 



where F = (ApB/ABB)(3ApBP^/2ep) is the interaction parameter of the system that embodies all interaction infor- 
mation necessary for understanding the zero-sound mode behavior. Alternatively, we can write 



F = 



/cpapB apB 
27r flB 



1 



mp 

TUB 



1 



mB 
mp 



(20) 



to expresses the interaction dependence in terms of the scattering lengths. 

We find that the angular dependence of the mode fluctuation is restricted by the density-density coupling of the 
liquids in the dilute fermion-BEC mixture. Expressing Wp in a spherical harmonics form, Up = lUimyim, we 
obtain that uim — for m ^ 0. In the collective excitation, the local Fermi momentum distribution fluctuates as a 
longitudinal mode with an azimuthal symmetry around the direction of propagation with the deformation amplitude 
Up, so that the sole momentum dependence is upon the angle between k and p. Rather than writing Up = uiqYiq, 
we simply determine the dependence upon p by solving for Up from Ea. H19|l and take the angular average. Defining 
the scaled phase velocity as s = uj/kvp, we obtain 



s^VF/cf- l + (fcO 



F 



s* 



-dx , 



(21) 
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where we set x = cos (p. From Eq. pil we can determine the dispersion relation of the collective modes. The complex 
conjugate, s* , follows from the appropriate choice of boundary conditions as we explain in the next section. 

The structure of the excited state is apparent from neither Eq. (|21|l nor Landau's semiclassical description (see 
[2^). We believe that the actual excited state can be written a linear superposition of two states: One is a product 
state of the ground state BEG and an excited fermion state, which is itself a linear superposition of particle-hole 
excitations as is the case in the pure fermi-liquid zero-sound mode. The other is the product of the ground state of 
the Fermi-liquid and a BEC-state that contains a single Bogoliubov quasi-particle excitation. 



IV. MODE ANALYSIS 



Equation 1)21(1 is a close relative of the one that describes the zero sound excitation of a two-component pure 
Fermi-liquid whose inter-particle interactions are restricted to s-wave partial waves |23| 

^--l + Slnfi±lV (22) 



Fo 2 \s-l 

where Fq is the Landau interaction parameter. 

When s > 1, the solution s to Eq. H21|l is real and the right-hand side of Eq. H21() reduces to a simple expression 



F 2 - 1, 

of the form of Eq. |(22J). The similarity suggests identifying F/ [s'^{vf/c)^ ~ + (^0^)] ^.s an effective Landau-liquid 
parameter. However, a true Landau-liquid parameter does not depend on the wavenumber k. On the other hand, the 
appearance of the wavenumber leads to the rich mode structure that we discuss in the following subsections. 

As mentioned before, F parameterizes all the dependence on the interaction parameters, and zero-sound modes 
exist for < F < 1. The eigenvalues of Eq. 122|l depend on the velocity ratio vp/c, which quantifies the ratio of the 
time scales of the boson and fermion dynamics. In fact, in the long wavelength limit, the oscillation periods of both 
the BEG and fermion modes of wavenumber k are order of {kc)~^ and (fcup)"^. These time-scales are also associated 
with the response time to density perturbations of spatial variation of k~^. 



A. Collective mode excitations < -F < 1 



If 3> c and < -F < 1, two collective modes exist: an undamped mode with a real s and a damped mode 
with a complex s. The imaginary part of s describes the decay rate of the damped mode and the zero-sound analogy 
reveals that the decay mechanism is indeed Landau damping: the creation of particle-hole pairs in the Fermi-system. 
Landau damping only occurs when s < 1 due to energy conservation: a particle-hole excitation with momenta K-|-k/2 
(particle) and K — k/2 (hole) near the Fermi-surface of momentum K has an oscillation frequency that is smaller 
than uj. In other words, KK. ■ \i/mF < vpk, so that a single particle hole excitation cannot absorb all the energy of 
the collective excitation. 



1. Undamped mode 

The scaled dispersion relations of the undamped modes, lo^/vf vs. fc^, at various velocity ratios of vf/c with a 
fixed F — 0.5 are shown in Figs. ^ 121 andO For vf > c the dispersion is linear with the phase velocity v « at 
small fc, similar to the long wavelength dispersion of the zero-sound mode in a weakly interacting pure Fermi-liquid. 
At larger k, the dispersion becomes identical to that of an elementary excitation of a single homogeneous BEG (i.e., 
Bogoliubove mode) with lo ~ kc\J\ + (fc^)^. The transition from linear dispersion to Bogoliubov dispersion occurs 
gradually around fctrans = ^^'^ \J {vp /cY — 1 as a result of the avoided crossing of the two dispersions (See Figs. ^ and 
121) On the other hand, when vf < c, the undamped mode dispersion is the Bogoliubov dispersion for the entire range 
of fc- values, as seen in Fig. |31 

Since our treatment of the Fermi- liquid dynamics assumes k <^ kp, we need to check the wavenumber regime 
of validity. However, for mixtures that consist of light boson and heavy fermion atoms (for example, ^H and ^"K, 
or ''Li and ^''Rb), this range can extend to wavenumbers comparable to or greater than as implied by fcp^ = 
(mp/2mB)('yp/c). 
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FIG. 1: The dispersion relation {lu£^/vf vs fc^) for a fixed F — 0.5 and a fixed vp/c = 2. The dispersion changes from linear to 
Bogoliubov dispersion. 
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FIG. 2: The dispersion relation {lu^/vf vs fc^) for a fixed F — 0.5 and a fixed vp/c = 3. The dispersion changes from linear to 
Bogoliubov dispersion. 

livp < c, the BEC responds sufficiently fast to follow the fermion oscillations adiabatically and the static description 
of BEC-mediated interaction is valid. In this case, the EEC-mediated interaction is attractive and there is no Fermi- 
liquid like zero-sound mode. Instead, there is a collective oscillation that resembles the elementary excitation of a 
pure BEC with v ~ c for all range of k (Fig. On the other hand, if vf ^ c, the BEC response is significantly 
slower than the fermion frequency and the BEC oscillates out of phase with the fermions. The overall sign of the left 
hand side in Eq. (|21l) is positive for small k, giving an effective repulsive fermion-fermion interaction. In this regime, 
the collective fermion mode with zero-sound characteristics occurs with phase velocity velocity v ~ v-p (see Figs. ^ 
and EI). 

2. Damped mode 

In addition to the undamped mode, a damped collective excitation exists in < F < 1. As stated before. Landau 
damping 5] is the dominant damping mechanism of zero-sound in the regime of collisionless dynamics. 

The imaginary part of the complex eigenvalues for Eq. (|21|l is proportional to the damping rate. Assigning the 
correct sign to the imaginary part of s in the eigenvalue equation Eq. f l21|l is rather delicate: While the damped mode 
corresponds io s — r ~ with r, 7 > 0, causality requires that the integral on the right-hand side, which describes 
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FIG. 3: The dispersion relation {lu^/vf vs fc^) for a fixed F — 0.5 and a fixed vp/c = 0.9. Tlie zero sound dispersion follows 
the Bogoliubov dispersion for all k values. 



the retarded fermion-mediated boson-boson interaction, acquires s* = s + ij. The damped excitation evolves in time 
as exp [—1 (w — iT/2)t] with T = 2kvFj. Independent from the collective mode equation, we estimate the damping 
rate r(A;) using the Fermi-golden rule. The initial state \i) is the product of the single quasi-particle BEC state with 
momentum k and a fermion state that consists of a completely filled Fermi sphere. The final state |/) is the product 
of the BEC-ground state and an excited fermion state with a particle-hole pair of momenta K -I- k/2 for the particle 
and K — k/2 for the hole. The quasi-particle-hole pair is created near the Fermi-surface of momentum K. Treating 
the fermion-boson density-density interaction as a perturbation Hp and summing over the final fermion states, we 
obtain the damping rate 



TMk) = i2n/h)Y,\{f\Hp\^)\'SiEf 
f 



Ei) = 



ttF 



kc 



(24) 



where the energy Ei and Ef represent the energies of the initial and final state, respectively. The damping rate is 
proportional to kc as well as F and (wp/c)"^. 

Note that Eq. (|24l) is valid for k <^ k-p. Figures 0] and [S] show Im[ijj]^/i;F obtained both by numerics from Eq. I|21() 
and by the Fermi-golden rule estimate for vp/c = 2, and F — 0.01 and F = 0.5, respectively. For < F ^ 1 with 
fp/c > 1, the rate of damping is orders of magnitude smaller than the excitation frequency, the phase velocity of the 
damped mode is almost equal to that of the BEC sound mode, and the Fermi-Golden rule approximation is in good 
agreement with the numerical solutions to Eq. H21|l . 

The damped mode exists for small wavenumbers k and ends at the 'termination point', fccnd, as seen in Figs.0]and 
|31 This behavior is reminiscent of the excitations in ^He ^2^. In the wavenumber interval that the mode exists, the 
rate for Landau-damping does not increase monotonically with k. This behavior is somewhat counter- intuitive as an 
increase in k generally increases the phase space volume of the available final states. The termination point /ccndC 
expressed in terms of 7 is 



fcond^ = (1 - 7') 



- 1 + 



F 

T 



7r7 



4 V2 



1/2 



(25) 



which can be tested experimentally. 

The magnitude of Im[cij] is significantly smaller than Re[cL;] even for a system of relatively strong interactions such 
as F = 0.5. Since we expect the collective excitations to be good quasi-particles, if the lifetime of the mode is short 
enough, the damping of the mode can be measured in atom trap experiments. The time scale of the mode damping 
is set by 



To = 



2vF 



mp 1 

TUB ckp 



(26) 



To is the time that takes a particle with velocity vp to travel half the distance of the BEC healing length. Alternatively, 

.01-1/3 ^ /'R_2^1/3|-^0ll/3 



we write fcp in terms of the average inter-fermion distance as rp = [pp\ 



(67r2 



L^fJ 



3.9/rp, and 
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FIG. 4: wf/c = 2, F = 0.01. The dashed hne shows the numerically evaluated Im[a;^/iiF], while the dotted line plots Im[Lj^/-i;F] 
evaluated from the Fermi-Golden Rule. 




0.5 1 1.5 



FIG. 5: vf /c — 2, F — 0.5. The solid line shows the numerically evaluated values of Re[tJ^/i;F], the dashed line plots the 
numerically evaluated curve of lm[uj^/vF], whereas the dotted line gives the values of Im[a;^/i;F] determined from the Fermi- 
Golden Rule calculation. 



find that tq is proportional to the time a particle requires to travel at the velocity of the BEC-sound in order to 
cover the average inter-fermion distance tq ~ (tof/S.Otob) x (r-p/c). Under typical experimental conditions tq can 
be a few milliseconds, which is sufficiently long to measure the damping and sufficiently short to keep a propagating 
wavepacket from leaving the trap region where the system is approximately homogeneous. 

Since the damped mode can be viewed as an elementary BEC-excitation that has been modified by the fermion- 
mediated boson-boson interactions, the decay can also occur by Belyaev damping where pairs of boson quasi-particles 
are created |2^. While our treatment does not describe this damping, the decay rate due to Belyaev damping tends 
to be slow in the long wavelength limit k < FbcI ~ (G^/tt/S) ^ P%(i%b i^O^^^^ thus can be neglected in the range 
of k of our interest. 

Finally, it is instructive to express s as a function of the velocity ratio, wf/c, in the limit fc ^ (see Fig. El) 
Similar to the behavior of the dispersion relations, there is a transition of the long wavelength phase velocity 
from its undamped mode v vp to the BEC-sound velocity v ~ c. The smoothness of this transition depends on the 
value of F: the smaller the value of F the more abrupt the transition occurs with varying velocity ratio. If v-p < c, the 
static description of the BEC-mediated fermion-fermion interaction is valid and the fermions experience an effective 
mutual attraction that does not support a long lived zero-sound excitation in a two-component pure Fermi-liquid 
system. We can determine the transition point (or region), {v-p / c)^^^.^^^^. For fc 0, the BEC and zero sound modes 
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FIG. 6: Plot showing the solution of Eq. 12111 for s at fc^ = as a function of cjvp ■ The solid line shows the undamped mode 
phase velocity with F = 0.01, the dash-dotted line plots its value for the undamped mode with F = 0.5, whereas the dashed 
line shows the values for the damped mode with F = 0.5 and the dotted line for the damped mode with F = 0.01. 

mix strongly near wp/c ~ 1, similar to the region where the linear dispersion merges into a Bogoliubov dispersion 
as we have seen in the dispersion relations. Quantifying the deviation of s from unity by setting s — I + e, where 
< e < 1 in the Eq.(EU we find 



c 











trans 





-1/2 



(27) 



The avoided crossing, termination point and strong mode mixing occur in the velocity ratio interval between (wp /c)^^.^^^ 
and 1. 

B. Instability and phase transition 

Equation (|21|l provides insight into the stability, or lack thereof, of a homogeneous fermion-boson mixture. If 
F > 1, s the eigenfrequency of the damped modes of longest wavelengths take on purely imaginary values. The purely 
imaginary value, and not just the acquisition of an imaginary part as in [2^ (a claim that was also corrected in |27]). 
signals the instability of the homogeneous mixture. Then Eq. f 121(1 is reduced to 



F 



1 + 7 arctan 



(28) 



Since the unstable mode is the lowest frequency mode, i.e., the damped or BEC-sound like mode in the stable 
mixture, we can say that it is the fermion-mediated boson-boson interaction that triggers triggers the instability of 
the homogeneous EEC in the mixture. Since F depends on fcp, the instability condition F > 1 can be expressed as a 
condition on the fermion density: p^p > pp^crit, where 



PF,ciit 



aBs/aFB 



(1 + mB/mp)(l -I- mp/me) 



(29) 



If the equilibrium fermion density exceeds pp^crit, the Fermi-liquid becomes immiscible to the EEC and spontaneously 
breaks the translational symmetry by undergo ing phase separation into cither (a) pure fermion and pure EEC or (b) 
pure fermion and mixed fermion- EEC phases 'll^. 

Equation (|28|) also provides insight into the dynamics of the onset (the early stages) of the phase separation 
instability. Consider an experiment that brings p p in R fermion-BEC mixturG to cibovG pF,crit 

by an abrupt change 

of obb if^ . In response, the amplitude of the longest wavelength modes of the damped zero-sound excitations grow 
exponentially at a rate 7(A;)fc^/ro nucleating spatial regions of pure fermion matter when the dynamics becomes 
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non-linear. Those modes that grow fastest (with wavenumber kd) dominate the dynamics and we expect them to 
determine the size of the single phase matter clusters and the rate of the phase separation p^. 

When 7 is sufficiently small, we can expand the right-hand side of Eq. 128|l and determine kd analytically. The 
modes with wavenumber in the range, € [0, are unstable and grow exponentially at the rate R — 'ykS^r^ 

with 



Rto = 

The maximum growth rate occurs at kd, where 

^ F - 1 ttF /wpN -2 




(l-F + q^) }k^ 



(30) 



128 V~/ 



(37rF)2 + 128 ( ^ ) (^^-1) 



(31) 



with the average cluster size ~ 27r/fcd. 



V. EXPERIMENTAL ASPECTS 



In this section, we mention relevant atom trap techniques and estimate magnitudes of quantities to gauge the 
feasibility of experimentally observing the collective mode physics. The important parameter for the boson-boson 
interactions in a dilute BEC is the dimensionless gas parameter, D — \/ Pb'^bb 122|- We find it useful to express the 
density in units of a reference density that corresponds to the lower range obtained in atom traps, pR — lO^^cm""^. 
Then, the gas parameter is 



D = ,/^^3.2xl0-^ (32) 
V PR Inm 

We introduce the atomic mass number of the bosonic and fermionic isotopes, and {Ay — 6 for ^Li, for instance) 
and express AbbPb ^^-'^ Fermi-energies e^? as 

6nK flBB 

AbbPb - — i "i 1 

Ab PR Inm 

.p=^r^V'\ (33) 



The time scale for the fermion dynamics is equal to 



2/3 

ip = — « 2.09psec x ^p ( ^ ) ., (34) 



from which it follows that the Fermi velocity can be of the order of several cm/sec, 

2 24.5 , /poy/' 
v-p — — w — — cm/sec x — . 35) 

Most importantly, the velocity ratio f p /c which quantifies the ratio of the relevant fermion to boson time scales and 
which determines the validity of the static approximation of the BEC-mcdiated interaction, can be written as 



VY _ Ab 0.55 ( p'^^ ^'^ 



c 



Ay i?i/3 



(36) 



Clearly experiments can access both the vp/c < 1 and up/c > 1 regimes. For example, ^Li with pp — lO-'^^cm"^ 
immersed in ^'^Na BEC with pg = 3 x 10"'^'*cm~'^ and cbb ~ 1.5nm gives vy/c « 0.46, whereas changing the density 
of the BEC to Pg = 3 X lO^^cm"^ increase the velocity ratio to vy/c~ 4.6. 

To observe the predicted collective modes an experimentalist could measure the group velocities in a fermion-boson 
mixture that is contained in a cigar shaped trap. We assume that the transverse confinement is not so tight as to 
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bring the system into the regime where it becomes effectively one-dimensional. A sudden perturbation near the trap 
middle created by a focused laser beam can excite both collective modes. In response, two wavepackets of different 
densities propagate outward with different group velocities. The density variations (fermion and boson) can be imaged 
directly, as in the first observation of the BEC-sound mode [s^- The slower wavepacket would correspond to the 
Landau damped mode with 

1 110 _i 
— « — — msec X 
To Af 




which can be slow enough to observe significant propagation while sufficiently fast to measure damping before the 
wavepacket reaches equilibrium densities that are significantly different from trap middle. 

In more sophisticated experiments, the initial perturbation can be controlled, and an analysis of the time evolution 
of the shape of the propagating wavepackets can be used to infer the dispersion relations and damping rates. In the 
past, the greatest control was achieved by accessing a two-photon resonance to cause low intensity Bragg scattering. 
In such an experiment the momentum of the excitation can be varied by changing the angle of the crossed laser 
beams. The laser beams need to be focused near the middle of the trap to probe the nearly uniform region of the 
fermion-boson mixture. Similar techniques have been used in 'Bogoliubov spectroscopy experiments' jsj]. 

Finally, we remark that by tuning the boson-boson scattering length close to its zero-point obb ~ 0, thereby 
lowers the critical fermion density can always trigger the phase separation phenomenon. By triggering the phase 
separation in a cigar shaped trap, the surface tension can prevent the single phase domains to move past each other, 
the experimenter can obtain a string of single phase 'droplets'. The size of these domains reflect the dominant mode 
wave vector. Such experimental study was demonstrated in the phase separation of multi-component BEC's [3^ . 

VI. CONCLUSIONS 

We have analyzed the collective modes of a single component fermion-BEC mixture in ultra-low temperature where 
the fermions behave as normal Fermi-liquid. 

If > c and < < 1, the homogeneous mixture is mechanically stable and two collective modes with zero-sound 
character exist. One is an undamped mode that resembles the zero-sound mode of a pure Fermi- liquid system. The 
dispersion of this mode is linear with phase velocity w/fc « up in the long wavelength limit and then merges into the 
Bogoliubov dispersion of an elementary BEC excitation near /ctrans = ^~^\/ (vf/c)^ ~ 1. The other mode resembles 
the sound mode of a pure BEC system in the long wavelength limit. The mode undergoes Landau damping and 
decays at the rate r(fc) = "f(k)kS,/To. The dispersion of this mode terminates at a well-defined wavenumber fccnd- 
In the long wavelength limit, the rate for Landau damping can be approximated using the Fermi-golden rule and is 
rFG(A:) = {7TF/2) {c/vF)ck. 

On the other hand, if < c and < F < 1, the BEC response is sufficiently fast to follow the fermion oscillations 
adiabatically, thus the dispersion is the Bogoliubov dispersion for all k. 

When F > 1, we have p^p > /Op^crit the homogenous mixture is unstable and undergoes phase separation by forming 
clusters of pure fermions. Our analysis shows that the growth rate is R = "fk^T^^ and the average cluster size is 
^ 2n/kd. 

Finally, our discussion of current atom trap technology and relevant observables indicate that the collective modes 
can be observed and the dispersions and damping rates may be measured. In addition, the phase separation can 
be triggered, enabling the phase separation dynamics characterization as soon as the fermions are cooled into the 
Fermi-liquid regime in the BEC-fermion mixture. 
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